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Abstract 

In this paper we study the energy spectra and related parallel momentum distributions of the neutrons from 
the breakup of 11 P>e. Earlier papers on transfer to the continuum reactions have shown that the breakup 
amplitude reduces to a simple eikonal approximation in the limit of weakly bound projectiles. The aim is to 
establish a reliable method to obtain information on the structure of the halo from the experimental results 
which depend on the reaction mechanism. Theoretical results for the widths of the parallel momentum 
distributions and for the total breakup cross section are compared to recent experimental data. 

PACS number(s):25.70.Hi, 21.10Gv,25.60Ge,25.70Mn,27.20+n 



1 



I. INTRODUCTION 



The fact that energy spectra of particles emitted to the continuum from a nucleus in 
a nuclear reaction bear some relation with their original momentum distribution has been 



known and studied for many years |]- [15| and the form and widths of the energy distribution 



have been related to the original momentum distribution of the emitted particle in the 
projectile. Reactions with final state having a neutron in the continuum have been studied 
in detail by us [[H]] and we arrived at the conclusion that there were dynamical effects 
modifying the original distributions | 20[j . An important consequence of our approach is that 
it is the component of the momentum distribution of the breakup neutron in the projectile 
parallel to the direction of relative motion that is responsible for the form of the energy 
distribution. This is in contrast to the results given in refs. P1,|2|,|TT| where it is the total 
momentum distribution of the neutron in the target which is important while it is similar to 
the approaches of |lO| , p!2|Jl3|| . Ideas and techniques very similar to ours have been recently 
used in [|1^] to the study of halo breakup. 

We will not repeat derivations given in earlier papers but, in order to make this paper 
more self contained, we will describe the physical motivation of our approach and the un- 
derlying assumptions. Our aim is to show that halo breakup can be studied by using simple 
analytical formulae which are more general and more accurate than the eikonal forms and 
reduce to them in the limit of zero initial binding energy. The main difference between the 
present approach and the eikonal approach comes from the fact that the rescattering of the 
breakup neutron on the target is calculated by the optical model. On the other hand the 
use of the semiclassical approximation for the relative motion of the two ions, discussed in 
the following, makes our method computationally much simpler than the DWBA type of 



approach used for example in Ref. [T5| 



One basic assumption is the masses of the target and projectile core are large compared 
with that of the neutron. Another is that the relative motion of the projectile and target can 
be treated classically. For the high energy reactions described in the present paper we assume 
that the center of the projectile core moves on a straight line path with constant velocity and 
with a particular impact parameter relative to the target. In our reference frame the target 
is at rest. We represent the interaction of the neutron with the projectile core by a single 
particle potential moving with the projectile velocity, and its interaction with the target by 
an optical potential. We use an approximate solution of the time dependent Schrodinger 
equation for the neutron moving under the influence of these two potentials with the initial 
condition that it is bound in a particular state in the projectile. The theory presented 
in Section II yields an expression Eq. (|2.2j ) for the momentum (or energy) distribution of 
the neutron in the final state as a function of the impact parameter. The cross-section is 
obtained by integrating over the impact parameter. 

One consequence of the use of a straight line trajectory with constant velocity for the 
motion of the projectile core relative to the target is that the theory does not satisfy the 
overall energy and center of mass momentum conservation conditions. The neutron receives 
energy from the time dependent potential fields associated with the relative motion of the 
target and projectile. For a given neutron final energy and momentum the overall energy and 
momentum conditions have to be put in by hand. We discuss this point at the beginning of 
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the next section. The use of straight line trajectories neglects Coulomb deflections. This is 
a good approximation for light targets and projectiles and high incident energies. We have 
discussed corrections for Coulomb deflections in earlier papers on transfer to bound 
states. The same methods could be used for breakup where necessary. 



II. THE BREAKUP CROSS SECTION 



The breakup theory presented in this paper evolved from a theory for neutron transfer to 
resonance states of the target ||16|| . In that case the final neutron energy Sf was the energy 
of the neutron resonance in the target. Thus in the present work £/ is the energy of the 
neutron relative to the target in the final state. For elastic breakup this is the same as the 
final laboratory energy of the neutron if the target recoil kinetic energy is neglected. In the 
case of compound nucleus formation e / is the excitation energy of the compound state above 
the neutron threshold in the residual nucleus. In the case of inelastic scattering it is the 
energy of the breakup neutron before it scatters from the target. This is equivalent to the 
sum of the excitation energy of the target final state and the final neutron energy relative 
to the target. If the target recoil kinetic energy is neglected the final kinetic energy Ef of 
the ejectile is given by the energy conservation condition 

Ef — E inc = Q = Ei — Ef, (2.1) 

where E inc is the initial incident energy of the projectile in the laboratory, Q is the reaction 
Q- value given by Q = Ei — Ef and £j is the initial neutron binding energy in the projectile. 
The projectile is supposed to remain in its ground state after the transfer. With this ap- 
proximation Eq. (|2.1| ) relates Ef to the energy loss spectrum of the ejectile. It is not possible 
to calculate corrections due to target recoil within the present theory because of the three- 
body character of the final state. This should be done to discuss experiments in which a full 
kinematical reconstruction is made. However such experiments are not available at present. 

The starting point of the present paper is the expression for the energy distribution of 
the breakup neutron given in Eq.(|2.2|) which was derived in ]IE|]. It is obtained by working 



in the target reference frame. The projectile-target relative motion is treated semiclassically 
by using a trajectory of the center of the projectile relative to the center of the target 
s(t) = d + vt with constant velocity v in the z-direction and impact parameter d in the 
xy-plane. 

Eq.( |2.2| ) gives the neutron transfer probability from a definite single particle state of 
energy £j, momentum j$ = y/—2m£i/h, and angular momentum /j in the projectile to a final 
continuum state of energy Ef, momentum kf = ^'IrriE f/h within an interval d,Ef. It is the 
sum of the transfer probabilities to each possible final //-state in the energy bin cfej 

dP_ 



X lf (\l-(S lf )\ 2 + l-\(S lf )\ 2 )B(l f ,k). (2.2) 



The physical interpretation is that the projectile brings up the neutron which is scattered 
into a continuum state by the target. The interaction of the neutron with the target is 
represented by the S-matrix for scattering of a free neutron by the target nucleus. 
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In the derivation of Eg. ( |2.2|) in [|16j the neutron final state was taken as a scattering 



wave function of energy ej, including all possible neutron-target final state interactions, and 
(Si f ) is the optical model, energy averaged S-matrix which describes the re-scattering of 



the neutron on the target [0. In Ref. |17J it was shown that the first term of Eq.(|23 
proportional to |1 — (Si f )\ 2 , gives the neutron elastic breakup spectrum while the second term 
proportional to the transmission coefficient T = 1 — \ (Si f )\ 2 gives the absorption spectrum. 
This term contains contributions from inelastic scattering of the breakup neutron by the 
target nucleus and also from compound nucleus formation. 

The factor B(lf, lj) is an elementary transfer probability which depends on the details of 
the initial and final states, on the energy of relative motion and on the distance of closest 
approach d between the two nuclei. Its explicit expression reads: 

where the arguments of the Legendre polynomials P^ and Pj. are respectively X^ = 1 + 
2(fci/7i) 2 and Xf = 2(k 2 /kf) 2 — 1. Also k\ = — (e, — Ef + \mv 2 )/{Tiv) and k 2 = — (e* — £/ — 
^mv 2 )/(hv) are the z-components of the neutron momentum in the initial and final state 
respectively, rj is the modulus of the transverse component of the neutron momentum. It is 
conserved during the breakup process and in fact it can be given both in terms of the initial 
or final neutron parameters as r] 2 = k\ + jf = k\ — k 2 . mv 2 /2 is the incident energy per 
nucleon at the distance of closest approach d for the ion-ion collision. 

In the case of a weakly bound projectile it is possible to obtain much simpler expressions 
for the elastic breakup and absorption probabilities. In fact k 2 which is the k z component 
of the neutron momentum with respect to target, in the limit of a very small initial binding 
energy takes the value k 2 ~ kf . Then the Legendre polynomial P;. ~ 1 since its argument 
is very close to one. By introducing the classical angular momentum A = k/b and 21 f + 1 = 
2(lf + 1/2) = 2A in Eq. ([2.2|) we can replace the sum over the final angular momenta with 
an integral over the neutron impact parameters b with respect to the target. Also when the 
neutron rescattering on the target takes place at relatively high energy (£/ — \vnv 2 \ the 
phase shift in Si. = e 2l5l s can be approximated by the eikonal form using a phenomenological 
optical potential such that Su ~ e~ %x ^ and x(b) = J^a ^(x, y, z')dz' where V 2 is a 
complex potential whose real and imaginary strengths are negative. Then 



dP(d) m 



de f h k 



oo 

bdb IC1 _ c -*x(&h|2 -L 1 _ \ e -ix(P)\2 







ft (d,h)\ 2 , (2.4) 



which is the product of the free neutron elastic (plus absorption) cross section, by the mod- 
ulus square of the initial state momentum distribution along the relative motion direction, 
for a fixed value of the distance d. The more accurate expression used in corresponds 
to replacing the d in the wave function ipi in the last term of Eq. ( p.4| ) by |b — d|. It can 
be obtained by approximating Pu by a Bessel function, Pi f (Xf) w I (2r)lf/kf) = Io(2rjb). 
The total breakup probability is the sum of an elastic plus absorptive term because the full 
final state wave function contains a S-matrix which is unitary. In instead we started 



directly from an eikonal wave function for the neutron final state and we obtained only the 
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elastic breakup term since in that case the S-matrix was not unitary. The results discussed 
in this paper were obtained with Eq. fl2.2Q , but Eq. (|2.4j ) is very useful to get an insight into 
the physics involved in halo breakup. 

In Eq.( |2.4j) \ipi\ 2 is the initial state momentum distribution [|2T] 



\Md,h)\- 



1 



21; 



-S m J^/ imi (d, h)\ 2 « |Ci| 



-2r\d 



2r]d 



(2.5) 



C\ is the asymptotic normalization constant of the initial wave function. The energy distri- 
bution of the final neutron has its main origin in the k\ dependence of ipi . 

The relation between the energy distribution probability and cross section is given by 



ClEf ■!() 



Pel(d) 



(2.6) 



as in Eq.(7) of P(sf,d) is given by Eq. fl2~2]) . P e i(d) = \Sct\ 2 is the ion-ion elastic 

scattering probability given by the modulus square of the projectile core-target S-matrix. 
The above factorized form of the neutron-target scattering by the core-target scattering has 



been used already in several papers on halo breakup , as Ref. |10],[T3],[l4],|2T| and and it is well 
known to hold for peripheral reactions [p3jp4f . In Ref. JTT] S ct was not factorised out of the 
breakup amplitude but since the results obtained both for the core momentum distribution 
and the neutron breakup cross section are consistent with those of Ref. [|l^] and with the 



results we present in this paper, the authors of [1 I ] conclude that the factorised form is a 



good approximation p5 |. 



The effect of \S ct \ 2 in the integral in Eq. Q2.6p is to cut-off small values of the core-target 



impact parameters. In Ref. |22| we used a parametrized form for \S ct \ 2 which corresponds 
to a smooth cut-off in d . The effect of the smooth cut-off depends on the parameter rj 
defined above. In the case of weak initial binding energy rj is very small and the smooth 
cut-off correction is negligible as compared to a sharp cut-off. Furthermore in a recent 



paper [15[ it has been shown that the breakup probability as a function of the core-target 
impact parameter, corresponding to P(ef,d)P e i(d) in our formalism, is peaked at a radius 
quite larger than the sum of the projectile-target radii and it rises sharply from zero to the 
maximum value (cf.Fig.2 of |15|). In |Lj| a DWBA approach has been used, which takes into 
account properly the distortion of the relative motion trajectory by the core-target optical 
potential. 

The above discussion justifies the use of a sharp cut-off in d. Thus a simple expression 
for the cross section is obtained, which is valid when P e \ is given by the strong absorption 
model such that P el = 1 if d > R x P p] = if d < R x [J21| 



el 

do _ c 2 S nRs 
dsj rj 



1 + 



dP{R s 
dEf 



(2.7) 



The above equation is consistent to leading order in - with the formulae given in [12 



for li = 0, 1. Eq.( [2.4|) is also similar to the breakup probability of [12] where instead of the 
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Glauber elastic plus absorption cross section factor the experimental neutron-target reaction 
cross section was used. The breakup cross section is sensitive to the choice of the strong 
absorption radius R s 



It can be estimated from 



1.4(4 /3 + ^ 3 )fm. 



(2.8) 



where A t and A pc are the mass numbers of the target and projectile core. 

The momentum distribution is one factor in the breakup probabilities expressions 
Eq. ( |2.2| ) and ( |2.4j ) but there are other factors which combine to modify its relation with the 
cross section Eq. (|2.7|) . One is the effect of the neutron-target rescattering shown explicitly 
by the Glauber factors of Eq.( |2.4| ). Another is the rj dependence in the factor ]L ^ L (l + ^r - ) 
in Eq.( [2.7|) coming from the ion-ion scattering and giving a kind of modified nucleus-nucleus 
geometrical cross section. 

The incident energy dependence in the total breakup probability is contained mainly 
in the neutron-target phase-shift Si f , as shown explicitly by the Glauber form x(b) ~ 
V 2 (b,0)/(vh). This is partly because of the factor 1/v in the expression for x(b) but also 
because the neutron-target optical potential V2 is energy dependent. However as we men- 
tioned after Eq.( |2.4|) the energy distribution of the final neutron has its main origin in the 
ki dependence of ipi(d,ki), the momentum distribution of the neutron in the initial state. 
With the approximate form Eq.( j2.5D it is contained mainly in the factor exp(— 2r]d). 



III. APPLICATION TO 11 BE BREAKUP 



The situation is especially striking for breakup of a halo nucleus like ll Be where the 
binding energy £j = —0.5MeV of the neutron in the initial nucleus is very small. The factor 
exp(— 2r]d) has a sharp maximum when 77 = 7^ = 0.155/m _1 or k\ = and Ef ~ mv 2 /2. The 
small value of 7^ gives distributions in k\ which are sharper than those obtained for normal 
heavy ions where the initial binding energy is of the order of —lOMeV and it makes the 
breakup cross section large. The values of rj as a function of the neutron final energy Ef are 
shown in Fig. (la) for the case of 11 Be by the solid line. Also shown is the same parameter for 
a 10 Be projectile (dot line) whose neutron separation energy is £j = — 6.8MeV and Zj = 1. 
There is a much stronger energy dependence and the minimum value of rj is much smaller 
in the case of 11 Be. The corresponding fci-dependence of ipi, calculated from Eqs. (|2.5|) is 
shown in Fig. (lb). As expected the momentum distribution is sharper in 11 Be (solid line) 
than in 10 Be (dashed line). It is also a characteristic of an s-state momentum distribution 
to be narrower than for p-states, d-states or larger Z-values. 

The initial state amplitude ipi(d, k\) depends on the choice of d = R s which in Fig. (lb) 
was taken to be 6.2 fm and which is appropriate for the breakup reaction of 11 Be on a 9 Be 
target. For a 208 Pb target we take R s = 11. 5 fm. These are close to values given by Eq. (|2.8| ). 
The widths at half maximum of the resulting momentum distributions are hAki = 49MeV/c 
for a 9 Be target and hAki — 39AMeV/c for a 208 Pb target. The difference is due to the 
change in the strong absorption radius. The dependence of the parallel momentum distri- 
bution on d or equivalently on R s can be understood by noting that the parallel momentum 
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distribution Eq.( [2.5|) is such that the square of its modulus represents the probability of 
finding a neutron in the projectile with a certain value of its momentum parallel to the 
relative motion velocity, when it is at a definite distance d from the center of the initial 
nucleus. 

The next step is to discuss whether information about the width of the neutron momen- 
tum distribution in the initial nucleus can be obtained from the breakup spectra. To be 
specific we focus on breakup reactions with a ll Be projectile on two targets. One possibility 
would be to use the measured spectra of breakup neutrons. In fact most of the present 
experiments use the energy loss spectra of the 19 Be fragments. If elastic breakup is the only 
mechanism leading to 10 Be fragments these two approaches would be equivalent. By en- 
ergy conservation, Eq. fl2.lp , the 10 Be energy loss spectra would have the same shape as the 
neutron spectra of Fig. (2). On the other hand the absorption component of Eqs.( |2.2| ) and 
( |2.4|) refers to neutrons which are detached from the projectile and are either inelastically 
scattered or cause some reactions in the target. They would have their energy degraded or 
might not been seen at all, but the energy loss spectrum of the ejectile can still be measured. 
In the figures the spectra are plotted as a function of Ef, but they can be interpreted as 
energy loss spectra using equation (|2.1| ). 

In Fig. (2) we show the energy spectra for the reaction 9 Be( n Be, 10 Be) 9 Be + n at 
Einc = 10, 41 and 72A.MeV calculated according to Eq.( [2.2| ) with an optical model S- 
matrix. The optical potential used is from Iffijj] . Technical details concerning this kind of 
calculations can be found in Refs. fT7jJl~8|j . For the 2si/2 initial state in 11 Be the param- 
eters appearing in Eqs. (|2.5| ) and (|2~7|) are: the spectroscopic factor C 2 S = 0.77, and the 
asymptotic normalization constant G\ = 0.94/m -1 / 2 . The dot lines in Fig. (2) are the elastic 
breakup spectra while the dashed lines are the absorption spectra. The solid lines give their 
sum which correspond to the ejectile inclusive energy spectra leading to 10 Be fragments. 
The relationship between the momentum scale in Fig. (lb) and the energy scale in Fig. (2) is 
contained in the definition of k\ given after Eq. fl2.3p . In the same way the neutron parallel 
momentum distribution after breakup can be obtained from Eq.( |2.7| ) by a simple change 
of variable from Ef to k\. Thus the simple expectation is that the width Aej of the neu- 
tron energy spectrum should be related to the width hAki of the halo neutron momentum 
distribution by 

hAki = Asf/v. (3.1) 



However because of the extra factors in Eq.( |2.7| ) there is a modification in the above 
relation. To show the extent of the modification Table I gives values of hAki obtained from 
the calculated energy spectra using the relation fl3.1p for a 9 Be target and a 208 Pb target. The 
optical potentials used were from [28[] and |29j respectively. In each column the first number 
refers to 208 Pb and the second to 9 Be . The numbers in the last column should be compared 
with the widths of the momentum distribution Eq. (|2.5|) , which is hAki — 39AMeV/c for 
the 208 Pb target, because R s = 11.5 and hAki = A9MeV/c for the 9 Be target from Fig. (lb), 
where R s = 6.2. In the case of Pb there is a close correspondence between the width of the 
calculated energy distribution and the width of the momentum distribution Aki. On the 
other hand the values of Aki extracted from the energy distribution for the Be target are 
20% smaller than the Aki from Fig. (lb). 
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To show this effect more clearly we have plotted in Fig. (3a) the differential cross section 
or energy spectra as a function of k\, at the same incident energies as Fig. (2). Diamonds, 
dot line and crosses are for E inc = 10 , 41, 72A.MeV respectively. The full line is the 
initial momentum distribution in 11 Be as in Fig. (lb). All curves are normalized at the 
same peak value. Fig. (3a) is for the Be target and one can see that the initial momentum 
distribution is wider than the values obtained from the energy spectra. For a Pb target 
instead, Fig. (3b) shows that the energy spectra at high energy reproduce well the initial 
momentum distribution width. The different scales in Figs. (3a) and (3b) reflect the difference 
in ipi(d, ki) when d changes. One possible reason for this behaviour is the fact that in the 
case of the Pb target the halo wave function is probed at a larger distance as compared 
to the Be case. Thus changing the target but keeping the incident energy fixed on can get 
information on the momentum distribution at different radii. 

Results which seem consistent with our discussion have been obtained by Kelley et al. 
pOfl who measured the parallel momentum distribution of 10 Be fragments on several targets. 
They found Aki = 41.6 ±2.1MeV/c on a Be target, in agreement with our value in Table I. 
Also their results show a 5% change in the widths depending on the target, again in agree- 
ment with the maximum variation shown by our calculations. A similar target dependence 
was remarked in [Rlfl in connection with 11 Li data. On the other hand the small Aki values 



obtained at low incident energies for the same targets are due to the fact that because of 
the relation between k% and Ef the tails of the momentum distributions in Fig. (lb) cannot 
be sampled by a low energy reaction. 

In the case of a heavy target the experimental spectra contain an important contribution 
from Coulomb breakup which we discuss in Here we anticipate that for the lead 

target the Coulomb breakup calculated spectra give HAki = 43MeV/c and 47MeV/c at 
E inc = HOMeV and 790MeV respectively. 

Widths calculated from the energy distributions obtained with the eikonal model Eq. ( |2.4| ) 
are about 10% smaller than the ones given in Table I obtained with an optical model final 
wave function using the more accurate energy distribution in Eq.( |2.2j ). The agreement 
between the two calculations improves increasing the incident energy and for heavy targets. 
However the absorption is always underestimated by the eikonal approximation. This is 
because the contribution from the very low partial waves which is correctly described by 
the optical modes S-matrix in Eq.( [2.2| ) is instead underestimated by the eikonal S-matrix. 
Elastic breakup is dominated by high partial waves which are well described by the eikonal 
S-matrix. In this respect it can be useful to remark that while for a light target partial 
waves up to If = 7 are enough to describe the neutron rescattering on the target, in the case 
of the Pb target we found that partial waves up to // = 25 have to be taken into account. 
This is possible and easy to do with the present method while it requires very cumbersome 
calculations in other methods like DWBA. 

The relative amount of nuclear vs. Coulomb breakup in the total breakup cross section, 
for heavy targets is currently a question of debate for the experimental and theoretical 



implications |33|| . In the case of heavy targets we are studying the important contribution 
to the exclusive cross sections from Coulomb breakup in The results could be used 

to disentangle the two mechanisms. Therefore we give in Table II the calculated cross 
sections integrated over energies and angles, for the breakup of ll Be on 9 Be, 48 Ti and 197 Au 
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34| at E inc = AlA.MeV and on a 208 Pb target |35| at E inc = 72A.MeV together with 



the experimental values for the inclusive and exclusive cross sections. The first row gives 
experimental and theoretical values for inclusive cross sections. The theoretical values are 
the sum of the elastic, Coulomb and inelastic contributions, given in the third to fifth rows. 
The second row gives the exclusive values. In this case the theoretical estimates are the 
sum of the Coulomb and nuclear elastic cross sections. The Coulomb breakup calculations 
were done in first order perturbation theory similarly to Our results agree well 



with experimental values and with other theoretical estimates [0,0 and show that for 
heavy targets the nuclear breakup, although smaller than the Coulomb breakup, cannot be 
considered negligible. The cross section for the lead target are slightly smaller than those 
for the gold target because the incident energy was higher. 

Our results for the reaction 9 Be( u Be, 10 Be) 9 Be + n were obtained with R s = 6.2/m. 
Using R s = 5.9/m gives a breakup = 0.416 while using R s = 6.5/m gives a breakup = 0.366. 
For comparison, the results for the neutron breakup from the l0 Be core at R s = 6.2/m 
are: a el = 0.0256 and a inel = 0.0376, where we assumed a spectroscopic factor C 2 S = 4 
for the lps/2 state in 10 Be. Since the experimental spectroscopic factor is in fact 2.1 and 
the rest of the occupation probability is spread over three excited bound states in 10 Be, the 
above value can be considered as a higher limit estimate for the core breakup contribution 
at Ei nc = HA.MeV. The elastic plus inelastic cross sections obtained using the eikonal 
formula (|2.4|) are smaller by about 25%. This is mainly due to the absorption term for 
which the eikonal approximation is not very accurate. Here and in the following we refer to 
the cross section coming from the absorptive term of Eq. (|2.2|) as inelastic. 

At this point it is worth mentioning that some authors |TT],|T2J refer to the absorptive 
term as to stripping, while others [13] have used the term inelastic breakup to include also 
processes in which the core nucleus interacts inelastically with the target, besides those 
discussed here which are due to the neutron-target interaction. However as it has been 
discussed in O] these processes are expected to give rise to a small core survival probability 
and therefore they have often been neglected. Since it is always the (A p — 1) nucleus which 
is detected, its possible excitations must be restricted to those below particle threshold. We 
have discussed and included such processes in our calculations for breakup of normal nuclei 
|i~8| , |20"H and the same can be done for halo projectiles. The values of cross sections for the 
core breakup given above give an estimate of such processes. Some experimental evidence 



for such reactions has recently been discussed by Hansen [37 



In view of future experiments in the incident energy range 3 — lOA.MeV it is useful to 
know the energy dependence of the nuclear elastic and inelastic total breakup cross section. 
In Fig. (4a) we show the total elastic (cross) and inelastic (star) cross sections and the sum 
(diamond) of them for a 9 Be target. They have been obtained from the integration of 
the energy spectra from Eqs.( [2.2|) and (|2.7|) calculated at different incident energy. The 
very interesting point is the relative behaviour at low energy of the two components of the 
breakup and also the dependence on the target. For this reason we show also in Figs. (4b) and 
(4c) the cross section values calculated for a 2S Si |38[ and a 208 Pb |35|] target. The neutron 
optical potentials used are from [^.^.^] respectively, for the three targets. Comparing 
the three cases we notice that at high incident energies the elastic and inelastic breakup 
tend to give the same cross section. This is the geometric limit of the nuclear cross section. 
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At low energies in the case of the light target 9 Be the elastic breakup is larger than the 
inelastic one and this pattern is maintained increasing the incident energy. This is because 
the neutron-light-target optical potentials have only a surface component of the imaginary 
part, due to the fact that these nuclei can have surface collective excitations but do not 
have enough density of levels to allow the absorption into compound nucleus. On the other 
hand in the case of 208 Pb the absorption is larger at all energies because there is a large 
probability of surface excitations leading to neutron inelastic breakup and also to the fact 
that the neutron-Pb optical potential has a volume term of almost constant magnitude at 
high energies. Our results are consistent with those of [|TT|,|i"2"[ i n the case of the Be target. 



IV. CONCLUSIONS 



To conclude, neutron breakup from the projectile in a heavy-ion reaction has been dis- 
cussed in this paper in the framework of the theory of transfer to the continuum reactions 
which uses an optical model final state wave function for the neutron and the relation with 
the eikonal model has been clarified. The latter is a good guideline to understand the mech- 
anism of halo breakup and it is quite accurate to describe the halo neutron breakup at high 
incident energy and on a heavy target. In other situations one should use Eq.( [2.2|) with 
an optical model S-matrix. Both methods can be used to calculate energy spectra but the 
angular distributions of the breakup neutron can be obtained only from the eikonal model 



| 21| which contains the explicit dependence on the neutron final momentum components, 



furthermore this method can be easily extended to relativistic energies, as it has been done 



in p 



The models have been applied to the breakup of 11 Be on several target nuclei. The 
neutron energy distribution has a peak centered at the incident energy per nucleon of the 
u Be projectile with a width related to the width of the momentum distribution of the halo 
neutron. To a first approximation there is a very simple relation between these two widths, 
but sample calculations show that, due to details of the reaction mechanism, there can be 
differences between the initial momentum distribution and the measured one which can 
reduce its original widths up to 20% for a light target nucleus. 

Integrated nuclear elastic and inelastic breakup cross sections have been calculated and 
compared to experimental values and to Coulomb breakup cross sections. In the case of 
heavy targets it is very important to estimate the two contributions separately. Our results 
agree well with experimental values and show that for heavy targets the nuclear breakup, 
although smaller than the Coulomb breakup, cannot be considered negligible. Finaly the 
incident energy dependence of the nuclear elastic, inelastic and total cross section has been 
studied. 

We are grateful to P.G.Hansen for comments on an early version of this paper and to G.F. 
Bertsch for discussions. One of us (A.B.) wishes to thank W.Tornow for discussions and for 
information on neutron-light-nucleus optical potentials and H.Horiuchi for the hospitality in 
Kyoto during the preparation of this paper. 
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Table I: Parallel momentum distribution widths for Be and Pb targets. The value of Hv is 
related to the incident energy but takes into account the effect of the nucleus-nucleus Coulomb 
barrier at R s 





As/ 


hv 


hAki 


(MeV) 


(MeV) 


(MeV • fm) 


(MeV/c) 


110 


4.4 5.0 


22.45 27.71 


38.6 35.5 


300 


8.8 9.0 


44.03 46.92 


39.3 37.8 


790 


15.0 15.5 


75.1 76.82 


39.3 39.8 


Target 


Pb Be 


Pb Be 


Pb Be 



Table II: Cross sections for Be, Ti, Au [M] and Pb [35] targets. The first row gives experimental 
and theoretical values for inclusive cross sections. The theoretical values are the sum of the elastic, 
Coulomb and inelastic contributions, given in the third to fifth rows. The second row gives the 
exclusive values. In this case the theoretical estimates are the sum of the Coulomb and nuclear 
elastic cross sections. 



a(b) 


Be 


Ti 


Au 


Pb 


( u Be, w Be) 
("Be, 10 Be + n) 
el 
Coul 
inel 


0.29 ±0.04 0.39b 
0.24 ±0.05 0.23 

0.22 
0.0074 

0.16 


0.65 ±09 0.70 
0.55 ±0.11 0.47 
0.28 
0.19 
0.234 


2.45 ±0.20 2.73 
2.5 ±0.5 2.37 
0.35 
2.02 
0.39 


- 2.39 
1.8 ±0.4 2.09 
0.3 
1.79 
0.3 



13 



Figure Captions 



Fig.l.(a) Values of the parameter rj from for ll Be solid line, and 10 Be dot line, (b) Initial 
momentum distribution for the 2si/ 2 state in n Be, solid line, and lp 1/ / 2 -state in 10 Be, dot 
line. 

Fig. 2. Final energy spectra from Eq.([E7|) for 9 Be( u Be, 10 BefBe + n at Ei nc = 10 , 41, 
72A.MeV . Dot lines are for elastic breakup, dashed lines for inelastic breakup and solid 
lines are their sum corresponding to the inclusive 10 Be spectrum via Eq.([2.1|). 

Fig. 3 (a). Final energy spectra as a function of k±, at the same incident energies as Fig. (2). 
Diamonds, dot line and crosses are for Ei nc = 10, 41, 72A.MeV respectively. Solid line is the 
initial momentum distribution as in Fig. (lb). All curves are normalized at the some peak 
value, (b) The same as Fig. (3a) for 208 P5( 11 Se, 10 Be) 20S Pb + n. 

Fig. 4 Incident energy dependence of the breakup cross sections on 9 Be (4a), 28 Si (4b) and 208 Pb 
(4c). Diamonds joined by full lines give the total breakup cross sections. Crosses with dot 
lines are for the elastic breakup while stars with dotdashed lines are the inelastic breakup 
cross sections. 
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s (MeV) 



-0.4 -0.2 0.2 0.4 -0.4 
k^fm 1 ) 



-0.2 0.2 0.4 

k^fm 1 ) 



